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ABSTRACT 

We use UV/optical and X-ray observations of 272 radio-quiet Type 1 AGNs and 
quasars to investigate the C iv Baldwin Effect (BEff). The UV/optical spectra are 
drawn from the Hubble Space Telescope, International Ultraviolet Explorer and Sloan 
Digital Sky Survey archives. The X-ray spectra are from the Chandra and XMM-Newton 
archives. We apply correlation and partial-correlation analyses to the equivalent widths, 
continuum monochromatic luminosities, and q ox , which characterizes the relative X-ray 
to UV brightness. The equivalent width of the C iv A1549 emission line is correlated 
with both a ox and luminosity. We find that by regressing Z„(2500 A) with EW(C iv) 
and a ox , we can obtain tigher correlations than by regressing ^(2500 A) with only 
EW(C iv). Both correlation and regression analyses imply that ^(2500 A) is not the 
only factor controlling the changes of EW(C iv); a ox (or, equivalently, the soft X-ray 
emission) plays a fundamental role in the formation and variation of C iv. Variability 
contributes at least 60% of the scatter of the EW(C iv)-Z„(2500 A) relation and at least 
75% of the scatter of the of the EW(C iv)-a ox relation. 

In our sample, narrow Fe Ka 6.4 keV emission lines are detected in 50 objects. 
Although narrow Fe Ka exhibits a BEff similar to that of C iv, its equivalent width 
has almost no dependence on either a ox or EW(C iv). This suggests that the majority 
of narrow Fe Ka emission is unlikely to be produced in the broad emission-line region. 
We do find suggestive correlations between the emission-line luminosities of C iv and 
Fe Ka, which could be potentially used to estimate the detectability of the Fe Ka line 
of quasars from rest-frame UV spectroscopic observations. 
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Introduction 



One of the important properties of AGNs is the relation between the emission-line strength, 
characterized by the equivalent width (EW) , and continuum l uminosity, becaus e it reveals that the 
regions emitting these two spectral components are associated. Baldwin! (|1977al ) found that the EW 
of C iv A1549 (C iv) is inversely correlated with the quasar monochromatic luminosit y at 1450 A, 



Z A (1450 A), namely, logEW(C iv) = fclogl„(1450 A) + b. ICarswell fc Smith! (jl97«l ) referred to 
this trend as the "B aldwin Effect" (B Eff), a designation now widely used to describe line strength- 
luminosity relations. Baldwin! (1 1977a! ) identified this relation using only 20 quasi-stellar objects with 
29.8 < log ^ (1450 A) < 32.0 and 1.24 < z < 3.53. Subsequent UV/optical surveys have enabled 
investigation of this relation with wider luminosity and redshift ranges (e.g., Kinney, Rivolo & Ko- 
ratkar 1990; Zamorani et al. 1992). It has been found that the BEff exists for no t only C IV but many 



other broad emission lines such as Lyq, C in] A1908, Si iv A1396, Mg n A2798 (iDietrich et al 



2001 



Van den Berk et al.l 12009. in preparation! ). UV iron emission lines (I Green. Forster fe Kuraszkiewicz 



2002 



and even forbidden lines such as [O n] A3727 and [Ne vl A3426 (|Croom et al J 120021 ) . Ap 



plying a spectral-composite technique (jVanden Berk et al 



2001 



) to the Sloan Digital Sky Sur - 



vey (SPSS; York e t al. 2000) Data Relea se Three (DR3) quasar catalog ([Schneider et al.ll2005l ). 
Vanden Berk et al.l (12009. in preparation! ) found that the BEff evolves with redshift, which is a 
source of scatter in this relation for a sample with a wide range of redshift. 

The X-ray BEff (or Fe Ka BEff), in which the EW of the narrow Fe Ka line at 6.4 keV 
(hereafter abbreviated to Fe Ka) is anti-correlated with X-ray luminos ity, l v (2 keV), was discovere d 
in the early 1990s from observations by the X-ray observatory Ginga ( 
This relation has be en subsequently confirmed using data fro m ASCA 



Nandra et al 



Jiang et al 



2006 



wasawa fe Taniguchilll 993') 



Tanaka. Inoue fe Hoi 



19971) and from Chand ra and XMM-Newton (jPage et alJbooi : lzhou fe Wang 



1994 



2005; 



Bianchi et al.ll2007l ). Possible sites of origin for narrow Fe Ka emission include 



the broad emission- line region (BELR), the outskirts of the accretion disk, and the molecular torus 
(e.g., Weaver et al. 1992, Antonucci 1993; Krolik, Madau & Zycki 1994). 

Although it is well accepted that the BEff exists for many UV/optical emission lines (e.g., 
Osmer & Shields 1999; Shields 2007), there is currently no theoretical model that provides a 
compelling and complete explanation of this well-known phenomenon. Several physical explanations 
have been proposed to account for the UV/optical BEff. 

One promising explanation is that the continuum shape may be luminosity dependent. In this 

model, the UV/optical BEff is due to the softening of the spectral energy distr ibution (SEP) at high 

luminosity, which lowers the ion populati ons having high ionization potentials (jNetzer. Laor fe Gondhalekar 



1992 



Korista. Baldwin fe Ferlandl Il998l ). It has been found, using E instein Observatory data 



that the quasar SEP, parameterized by a Q Jl| ( Tananbaum et al. 1979 ). depends on UV luminos- 



1 Defined as a ox = log [l v {2 keV)/Z„(2500 A)] / log [v{2 keV)/i/(2500 A)] = 0.3838 log [l v {2 keV)/J„(250Q A)]. a ox 
is used to characterize the spectral hardness in the UV to X-ray band (e.g., Avni & Tananbaum 1982, 1986; Anderson 
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ity (e.g., Zamorani et al. 1981). Later resea rch using radio-quiet ( RQ) optically selected quasar 
samples from the SDSS Early Data Release (jStoughton et al.ll2002l ) confirmed and extended this 
result (e.g., Vignali, Brandt Schneider 2003). Using optically selected AGNs IStrateva et al. 
(|2005l ) and ISteffen et al.l (|2006l ) firmly established the correlation of a ox with UV luminosity for 



these sources. The idea that the UV/o ptical BEff is attribut a ble to SED-driven ionization ef- 



fects is supported both observationally (jZheng Malkan 



1993 



Greenl 11999 . Il998h and theoreti- 



cally (jNetzer. Laor fe Gondhalekarlll992l ). Recent work on a sample of non-Broad Absorption Line 
(BAL), radio-quiet, optically selected quasars indicates that the EW of C iv depends both on UV 
and X-ray luminosity. The physics of the C iv BEff is apparently associated with both UV and 
X-ray emission (e.g., Gibson, Brandt & Schneider, 2008, hereafter GBS08, and references therein). 

Other proposed BEff drivers include the Eddington ratio , L/L-^dd ( Baskin fe Laorll2004 : 



2004 



1992 



Warner. Hamann fc Dietrich 



Wandel . Peterson fe Ma lkan 



ity (jWarner. Hamann fe Dietrich! 120041 ) . 



2004 



1999 



Shields 



s 2007), 



Bachev et al 



Zhou fc Wanlh oOfih. the black-hole mass (jNetzer. Laor fe Gondhalekarl 



and the luminosity dependence of metallic- 



In this paper we investigate the origin of the BEff for the C iv emission line in a sample of 
272 Type 1 AGNs and quasars. Although C iv is not the only UV/optical broad emission line that 
exhibits a BEff, we selected it to study this phenomenon not only because it is a representative 
and well-accepted BEff emission line, but also because C iv resides in a relatively clean spectral 
region where the local continuum can be well approximated as a single power-law, with few blends 
with other emission lines (in particular the iron emission forest) and limited contamination from 
the AGN host galaxy. These properties make it relatively straightforward to perform spectral 
fitting and obtain accurate emission-line parameters for C iv. We also use partial-correlation 
analysis (PCA) and linear-regression regression analysis to investigate the correlations between 
EW, monochromatic luminosity, and a ox for C iv and narrow Fe Ko emission lines. 

Over the past three decades, there have been a large number of studies of the BEff. Our work 
on the C iv/Fe Ka lines combines the following important features (1) a wide range in redshift 
(0.009 < z < 4.720) and luminosity (27.81 < logZ„(2500 A) < 33.04) that allows one to disentangle 
evolutionary vs. luminosity, so that we are not narrowing our study for quasars with a particular 
luminosity or at a certain redshift; (2) a relatively large sample size (272 objects); (3) the use of 
partial correlation analysis; and (4) a high X-ray detection rate (~ 94%); and (5) estimates of the 
effects of observational errors and object variability. 

We describe the sample selection in § [2] and the methods used to process the data in § [3j In 
§ HI we perform partial-correlation and linear-regression analyses to investigate the roles of a ox in 
the C iv and Fe Ka BEffs. In § [5j we probe the connections between C IV and Fe Ka relationships 
in EWs, fluxes and luminosities. We present our conclusions in § [6J Throughout this work, we 



& Margon 1987; Wilkes et al. 1994; Vignali, Brandt & Schneider 2003; Strateva et al. 2005; Steffen et al. 2006; Just 
et al. 2007). 
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adopt the following cosmology: = 0.3, Qa = 0.7, Hq = 70 km s Mpc 



2. Sample Construction 



The quasar sample in our study is drawn from three sources: 50 objects from Ijiang et al 



(J2006|), which will be referred to as "S ample A"; 9 8 obi ects from GBS08, which will be referred to 
as "Sample B"; and 124 objects from lJust et al.1 (|2007l ). which is referred to as "Sample C". We 



define our "combined sample" as the combination of these three data sets. 



2.1. Sample A 



Jiang et al.l (|2006l ) compiled a dataset of 101 Type 1 AGNs with Fe Ka observations from both 
the Chandra and XMM-Newton archives. The detection fraction of the Fe Ka line in their combined 
sample is around 55%. The redshifts of the AGNs range from 0.003 to 3.366, but most of the objects 
(87%) are low-redshift (z < 0.4) AGNs. The monochromatic luminosities, ^(2500 A), range from 
10 26 ' to 10 31 ' 5 erg s^ 1 Hz -1 . We chose this dataset because it is the most complete Fe Ka BEff 
sample with high-quality data obtained from the most sensitive X-ray missions. However, because 
the number of AGNs observed in the X-ray band is much smaller than the number observed optically 
and not all X-ray observed AGNs present Fe Ka emission lines, the Fe Ka sample is significantly 
limited in size. 

We searched for UV/optical spectra of all 101 objects from the archival databases for HST and 
IUE. We found 82 spectra covering the wavelength region around the C iv emission line (containing 
at least the 1500-1600 A band). If observations are available from both HST and IUE, we selected 
the HST observations due to their generally higher signal-to-noise (S/N) ratio. For spectra with 
multiple observations using the same instrument, we preferentially use spectra with higher S/N. 

Next, we excluded all the radio-loud (RL) objects from the core sample, because additional 
X-ray emission is produced by the radio jet (e.g., Brinkmann et al. 2000) and changes the value of 
a ox as well as the slope of the Fe Ka BEff. 

We further excluded 3 objects for the following reasons: 

• MCG-06-30-15: The Fe Ka profile of MCG-06-30-15 is well fit using a broad disk line model 
(e.g., Tanaka et al. 1995). The narrow component is not well resolved or very weak. In this 
work, we only study the narrow component of Fe Ka, so this object is excluded. 

• IC 4329a: The spectrum has low S/N, and the C iv emission line is cannot be accurately 
measured (e.g., Crenshaw & Kraemer 2001). 



PG 1407+265: This object was termed as an "unusual" quasar ([McDowell et al.lll995l ) be- 



cause it contains extremely weak Lya and C iv lines. Because of its peculiarity, we exclude 



- 5 - 



it from our sample. 



We removed 5 objects wit h strong associated ab sorption li nes of C IV A 1549: P G 1411+442 



( Wise et alJl2004h. N GC 4151 (ICrenshaw et al.lll999lh Ark 5 64 (|Crenshaw et al.lll999T ). NGC 4051 



Collinge et al.ll200ll ). and PG 1114+445 (jShang et al.1 120071 ). These features prohibit reconstruc- 



tion of the unabsorbed C IV A1549 emission-line pro file. The X-ray absorption associated with 
the UV line absorption (jBrandt. Laor 8z Wills! |2000| ) might also lead to an under-prediction of 
continuum flux at 2 keV and, therefore, an incorrect estimation of the intrinsic value of a ox . 

Finally , we exclude 8 objects that are classified as Seyfert 1.5 (Sy 1.5) and Sy 1.90 ( Osterbrock 



1981 



19891 ) . These objects are intermediate between Sy 1 and Sy 2 galaxies and are often subjected 
to obscuration along the line of sight. Thus their UV, X-ray luminosities and a ox values are also 
potentially affected. 

The final version of Sample A consists of 50 AGNs. Among these objects, 34 are found in the 
HST archive (FO£@ or STISj), and 16 are found in the IUE archive (SWF@ or LWF0). The Fe Ka 



detection rate is 55%, the same as in the entire Ijiang et al. study. 



2.2. Sample B 

Objects in Sample B were selected from 536 SDSS DR5 (Adelm an-McCarthy et al. 2 007) 



quasars in GBS08. These quasars, taken from the DR5 Quasar Catalog ([Schneider et al.ll2007l ). are 
at redshift 1.7 < z < 2.7 and have been observed by Chandra or XMM- Newton. The lower redshift 
limit ensures that all SDSS spectra in this sample cover the C IV region; the upper redshift limit 
ensures that the rest-frame flux at 2500 A is covered so that a ox can be measured accurately. 

We excluded the BAL quasars in Table 1 of GBS08. BAL quasars can have strong absorption 
features in the C iv spectral region that prohibit accurate fitting of the continuum and emission- 
line profiles. In addition, BALs are usually associated with relatively strong X-ray absorption 
(e.g., Brandt, Laor & Wills 2000; Gallagher et al. 2006). We only retain objects with Chandra 
observations with angular offsets < 10' to avoid large X-ray flux uncertainties caused by variations 
of the point spread function. This restriction reduces our sample size to 149. In addition, we 
excluded RL objects and strong associated absorption-line (AAL) objects. The final version of 
Sample B consists of 98 objects with ^(2500 A) between 10 30 53 and 10 3L67 erg s _1 Hz -1 . This is 



2 Based on NASA/IPAC Extragalactic Database: http://nedwww.ipac.caltech.edu/ 

3 Faint Object Spectrograph 

4 Space Telescope Imaging Spectrograph 

5 Short Wavelength Prime 



(i 



Long Wavelength Prime 
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a relatively narrow luminosity range; however, it does provide a significant extension of Sample A 
since the latter is mostly composed of low-luminosity and low-redshift AGNs. 



2.3. Sample C 



To examine the relations between a ox , ^(2500 A), and l u {2 keVl. l.Tust etaJ d2007h com piled 
a sampl e of 372 objects , inclu ding 26 from their core sample, 332 from IstefFen et ah ( 2006 ). and 
14 from lShemmer et al.l (|2006l ). BAL quasars, RL quasars, and gravitationally lensed objects have 
already been excluded from this sample. The gravitationally lensed quasars are removed because 
their fluxes are strongly amplified and thus their luminosities are uncertain. 

Among the 372 objects, 38 are already in Sample A. For the rest of the AGNs, we searched 
for existing spectra with C iv coverage preferentially from SDSS, then the HST and IUE archives. 
Finally, we removed 5 AGNs whose spectra contain strong AALs. These restrictions leave 124 
objects in Sample C, in which 91 objects are from the SDSS DR5 quasar catalog, 13 from the HST 
archive, and 20 from the IUE archive. The redshift of this sample ranges from 0.015 to 4.720 and 
Z„(2500 A) ranges from 10 28 - 12 to 10 32 32 erg s" 1 Hz" 1 . 



2.4. Combined Sample 

The combined sample (Table [TJ consists of a total of 272 objects: 189 (69.5%) have spectra 
from SDSS, 47 (17.3%) from HST, and 36 (13.2%) from IUE. The redshifts range from 0.009 to 
4.720 (Fig. [TJ). The gap between z ~ 0.5 and z ~ 1.5 is caused by instrumentation limitations. 
Because of the wavelength coverage of the SDSS spectrographs, the redshifts of SDSS quasars 
having C iv coverage must be greater than 1.5. Most intermediate-redshift AGNs (0.5 < z < 1.5) 
are too faint for their UV/optical spectra to be taken by IUE and HST. Our sample exhibits a 
strong redshift-luminosity correlation (Fig. [1]); we discuss this issue further in our analyses below. 

The Z v (2500 A) of this combined sample ranges between 10 26 - 53 and 10 33 - 04 erg s" 1 Hz" 1 , 
including Seyfert galaxies to the most-luminous quasars in the Universe. The X-ray detection rate 
is 94.9%. The UV properties of the combined sample as well as the UV and X-ray properties of 
Sample A are tabulated in Table [2] and Table The combined sample will be used to investigate 
the C iv BEff, while only Sample A will be used to study the Fe Ka BEff. 
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3. Data Processing 
3.1. Bad Pixel Removal and Reddening Correction 



To ensure that we use high-quality data to perform the continuum and emission-line fitting, 
we remove bad pixels in the SDSS spectra based on the mask column contained in the SDSS quasar 
spectral files. We removed all the bad pixels in the spectra of SDSS objects in Sample C|H The 
excluded pixels cover less than 10% of the total pixels for over 98% of SDSS objects, and the 
maximum fraction of removed pixels for a single object is 15%. 

We perfor m Galactic reddenin g corrections to all the spectra using the E( B — V) dependent ex- 



tinctio n curve of lFitzpatrick! (|1999l ). Values of E(B— V) are calculated following lSchlegel. Finkbeiner k, Davis 
(|l998h . 



3.2. Spectral Fitting 

3.2.1. Sample A 

For each UV/optical spectrum from HST or IUE, we fit the local continuum in the vicinity 
of C iv (typically 1300-1700 A) using a single power-law. We do not expect our measurements to 
be significantly affected by host-galaxy components because 1) our sample contains only Type I 
AGNs and quasars in which emission from the nucleus dominates the light from host galaxies in 
the UV band; and 2) C iv is in the UV region in which the non-nuclear emission primarily arises 
from massive stars such as O and B stars. Examination of the spectra does not reveal any stellar 
absorption lines, indicating that the contribution from starlight is negligible. We do not subtract 
the iron emission forest as this component is usually not strong around the C iv emission line (e.g., 
Shen et al. 2008), and the wavelength coverage of the HST and IUE spectra is frequently too 
narrow (only a few hundred Angstroms) to fit this component. The Balmer continuum "small blue 
bump" only appears in the wavelength range between 2000-4000 A, so its contribution is negligible 
around the C iv region. 

The emission-line spectrum is obtained after subtracting the power-law continuum. We fit the 
C iv emission lines using two Gaussian profiles. The model always produces visually acceptable 
fits. We mask narrow absorption-line features appearing near the emission-lines so as not to under- 
predict the emission line flux. The equivalent widths, the emission-line luminosities under the 
assumption of isotropy, and the continuum monochromatic luminosities at 2500 A are calculated 
under our adopted cosmology and are tabulated in Table [2j 



To first order, we use 2-10 keV luminosities tabulated in Table 1 of iJiang et al.l (|2006l ) to 



7 For data processing of Sample B, refer to GBS08. 
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estimate l v {2 keV) under the assumption that the X-ray continuum is a single power- law with 
photon index T = 2 from 2 keV to 10 keV (e.g., Page et al. 2005; Shemmer et al. 2005; Vignali et 
al. 2005), so that 

L(2-10 keV) 



l u (2 keV) 



(1) 



i>2 In 5 

where hv<i = 2 keV and h is Planck's constant. To obtain the Fe Ka emission-line flux, we further 
assume that the Fe Ka emission line resembles a single Gaussian profile: f\(v) = A ■ e~( v ~ uo ^ l 2a ' 2 . 
This allows one to express the line flux F\ in terms of l u {2 keV) or ^(2-10 keV), the Fe Ka 
equivalent width EW(Fe Ka), and the central energy (eo = 6.4 keV) of the emission line. Assuming 
the continuum flux f c = C ■ u^ 1 , we can derive 



EW(Fe Ka) 



dv 



The Fe Ka emission-line flux is 



A 

C 



2na 



- x2 ' 2 ° 2 dx 



'A) 



A 
C 



Ft 



h(y)dv = A 



r+co \ 

V2^a - / e ~ x2/2a2 dx w AV^a. 

J l/Q J 



2ttvqo. 
(2) 



(3) 



The approximations are valid because generally vq S> a so Vq/2ct 2 3> 1. For instance, eo(Fe Ka) = 
6.4 keV while the width, <r(Fe Ka), is usually < 0.1 keV; result, e~ u o/^ 2 « 0. We can therefore 
calculate the Fe Ka line flux by 



F, 



EW(Fe Ka) F(2-10 keV) 



In 5 



(4) 



The Fe Ka EWs, emission-line luminosities, and 2 keV monochromatic luminosities of Sample A 
are tabulated in Table O 



3.2.2. Sample B 

For objects in Sample B, we directly adopt the fitting results from GBS08. In their paper, the 
SDSS spectral continua were fit with polynomials, and the C iv emission lines were fit with Voigt 
profiles. The different model in GBS08 from our work used to fit the C iv spectral region will not 
cause significant differences; because the continuum around C IV is not contaminated with other 
emission/absorption lines, the polynomial fit will produce almost the same result as the simple 
power-law fit. In addition, since both multiple Gaussian and Voigt profiles produce acceptable fits 
to the emission line, they will give nearly the same line flux. The X-ray spectral continua were fit 
using a broken power-law with the power-law break fixed at 2 keV in the rest-frame in order to 
obtain l v (2 keV). 
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3.2.3. Sample C 



We fit the UV/op tical spectra from the SDSS using a routine, described bv lVanden Berk et al 
(|2009. in preparation! ) . developed for SDSS quasar spectra. This routine fits the "underlying contin- 



uum" using three components simultaneously: a power-law, the iron emissio n forest, and the small 
blue bump. We adopt the UV iron emission templat e (1075-3090 A) from IVestergaard &: Wilkes 
feoOll ) and the optical template (3535-7534 A) from Iveron-Cettv. Jolv fc Veronl fcoO-d ). We first 
make a preliminary estimate of the power-law component by connecting two "line-free" points in 
the spectra. This provides initial estimates of the power-law parameters for the subsequent com- 
prehensive processing in which the spectra are fit by considering all three components mentioned 
above. We evaluate the fitting quality by calculating \ 2 values within some "line-free" windows. 
Finally, the continuum fit is subtracted from the original spectrum and the residuals are used to 
conduct emission-line fitting. For Sample C, the C iv emission-lines are fit by superpositions of 
two Gaussian profiles, which always yields acceptable fits for the data. We then calculate the C iv 
EW, emission-line luminosity, and monochromatic lumi nosity at 2500 A ( Table [2]). We adopt the 
values of a ox an d 2 keV monochromatic luminosity from bust et al.1 ( 2007 ). 



Examples of continuum and emission-line fits of the HST, IUE and SDSS spectra are presented 
in Fig. [2J 



3.3. Error Analysis 

In order to estimate the uncertainties in the measured quantities for objects in Samples A 
and C, we ran Monte Carlo simulations assuming a model with a perfect correlation between 
EW(C iv), /(C iv), and 7^(2500 A). For each spectrum, we add random noise to the original 
best fit to produce artificial spectra. The random noise follows a Gaussian distribution, and its 
amplitude is determined in one of two ways. For a spectrum from the HST or IUE database, we 
apply a low band pass filter, filtering out low-frequency signals via fast Fourier transformation! 
The residual signal is mostly noise. We then calculate the root-mean-square (RMS) of the noise 
and take this value as the random noise amplitude to be added onto the model spectrum. For a 
spectrum from the SDSS database, we simply use the uncertainty level associated with each pixel 
as the noise amplitude. For each spectrum, we produce 100 artificial spectra and fit them in exactly 
the same way as the observed spectrum. The error bar of a spectral parameter is then calculated 
as the RMS of 100 fitting results. The typical error bar shown in each plot (e.g., Fig. [3]) is the 
median of all the error bars of points in that plot. 

When evaluating the uncertainties of monochromatic luminosities, e.g., ^(2500 A), we must 
consider the contribution from the error in the luminosity distances (important for the low-redshift 



http: / /www.msi. umn.edu / software / idl/tutorial / idl-signal.html 
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objects). To validate this, we calculate the ratio of luminosity uncertainties without considering 
distance errors (Sf, the error in the flux measurement) to luminosity uncertainties after considering 
distance errors (5/,d) for objects in Sample A, and denote it as <5//<5/,d, m which / stands for flux 
and d stands for distance. Assuming that the square of total uncertainty can be expressed as the 
quadratic summation of the uncertainties of flux and distance individually, the square of the ratio, 
(Sf/5f : d) 2 , is more relevant than the ratio itself. We find that there are 5 out of 50 objects (in 
Sample A) whose ratios are above 0.1; three of them are even greater than 0.5. These objects have 
very low redshifts but large redshift uncertainties (5 Z > 0.001). For consistency, it is necessary 
to include the distance into the gross uncertainty calculation for all the objects in our samples. 
These uncertainties are also estimated using a Monte Carlo method. We adopt the most accurate 
values of redshift and their uncertainties from the NASA/IPAC Extragalactic Database^ These 
uncertainties are treated as noise amplitudes to be added to the redshift values. The luminosity 
uncertainties are calculated using the following error-propagation equation: 

SL = i7rd L ^(d L 5f) 2 +4/2 (5d L ) 2 

We adopt a 20% uncertainty for each X-ray continuum luminosity, e.g., l u (2 keV); their mea- 
surement errors are not available in the literature. The relative uncertainty varies considerably 
depending on the number of X-ray counts. For Sample A, when the Fe Ka emission lines are 
detected, the total number of counts is at least ~ 500. The Fe Ka emission-line luminosity and 
flux errors are calculated using the maximum error estimates. For instance, the upper bound of 
L(Fe Ka) is calculated using the upper bounds of both EW(Fe Ka) and L(2-10 keV); the upper 
bound of /(Fe Ka) is calculated using the upper bound of L(Fe Ka) and the lower bound of lumi- 
nosity distance d\, — 5d^. This uncertainty ignores any systematic uncertainty produced by errors 
in the cosmological model. 

For Sample B, because the parameter uncertainties are not given in GBS08, we simply apply 
a 20% uncertainty for all luminosity values as a first-order approximation. 

4. Drivers of the C iv Baldwin Effect 

4.1. Comparison with Previous Work 

We will first examine some important relations to determine if our measurements are consistent 
with previous work. It has been argued that a ox has no detectable redshift dependence (e.g., 
Strateva et al. 2005; Steffen et al. 2006; but see Kelly et al. 2007), so in this paper we neglect any 
redshift evolution of a ox - 



' ht t p : / / ned www .ipac.caltech.edu/ 
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Fig. [3] displays the plot of EW(C iv) against ^(2500 A) for the combined sample, distinguished 
by luminosity. We use the monochromatic luminosity at 2500 A rather than the traditional BEff 
wavelength (1450 A) because our choice is more convenient to compare the BEff with the correlation 
between C IV and a ox - The luminosities at these two wavelengths are well correlated. The quantity 
/1400//2500 is Gaussian distributed with a dispersion of a ~ 0.15 (Fig. 3 in Gibson et al. 2009), so 
using / I/ (2500 A) instead of ^(1450 A) should only add a small dispersion to the data points but will 
not significantly affect the slope of the BEff. We fit the data points linearly in logarithmic space 
using the EM (Expectation-Maximization) method (Dempster, Laird & Rubin 1977; Table H}. It 
is clear that EW(C iv) decreases with ^(2500 A) (Fig. 

It has b een reported that the slope of the BEff becomes steeper for high-luminosity quasars. 
For example, iDietrich et al.l (|2002l ) obtained a C iv BEff slope (—0.14 ± 0.02), which was shallower 
than the value reported in previous studies (—0.22 ± 0.05, Green 1996; also see Osmer, Porter & 
Green 1994; Laor et al. 1995). Using only t he EW(C iv) measur ements for their high- luminosity 



subsample with AL^(1450 A) > 10 44 erg s _1 , IDietrich et al.1 ( 20021 ) obtained a steeper slope of the 
BEff of —0.20 ± 0.03 for C iv. To investigate the slope change with luminosity, we fit our data 
points with /^(2500 A)< 30.5. We find that the slope of the low-luminosity sample is consistent 
with the slope of the entire sample wit hin 1 a (Table II). Th erefore, although our dataset exhibits 



(|2002l ). we do not find significant changes of 



a suggestive trend which disagrees with IDietrich et al. 
slope over luminosity. 

The large scatter in the BEff could have several causes, including observational error, intrinsic 
variation of the BEff (Osmer &; Shields 1999; Shields 2007 and references therein), luminosity 
dependence of the BEff slope, and the redshift dependence of the BEff. Vanden Berk et al. (2009, 
in preparation) found that the slope of the BEff does not change across redshift, but its scaling 
factor (or equivalently, the EW of a broad emission line at a fixed monochromatic luminosity) 
exhibits a significant change Our BEff slope is an overall average across the redshift and luminosity 
range we covered; remember that there is a strong redshift-luminosity correlation in our sample. 

We examined the ^(2500 A)-Z„(2 keV) and Z„(2500 A) -a ox relations for both Sample A and the 
combined sample using survival analysis (ASURV, Lavalley, Isobe & Feigelson, 1992) if censored 
data are involved, and find that all the results are consistent with previous work. When we calculate 
the slope of the Fe Ka BEff, we apply the Buckley- James method (Tabled Buckley <fc James 1979; 
Lavalley, Isobe & Feigelson 1992) included in ASURV to perform linear regression because the 
EW(Fe Ka) contains censored values. The slope of the Fe Ka BEff o f Sample A ( —0.11 5 ± 0.062) 
is consistent with the result from the RQ sample (—0.102 ± 0.052) in I Jiang et al.l (|2006l ). 
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4.2. The Effects of a ox on the C iv and Fe Ka BEffs 

Fig. U] shows the plot of EW(C iv) against a ox for the combined sample; the regression result 
from the EM algorithm for the linear relation is 

log EW(C iv) = (1.035 ± 0.075)a ox + (3.301 ± 0.119) (5) 

and the Spearman correlation coefficient is 0.607 (P < 0.001)0 Because a ox is an indicator of 
the hardness of the SED which controls the ionization level of C iv surrounding the central engine, 
Fig. H] demonstrates that as the ionizing flux becomes harder (a ox increases), the C iv emission has 
a strong positive response to a ox . 

Both a ox and Z;,(2500 A) are correlated with EW(C IV); which is a more fundamental driver? 
To investigate this issue, we applied PCA to EW(C iv), ^(2500 A), and a ox using the combined 
sample, Sample A, and a reduced sample. Table [5] presents Pearson, Spearman, and Kendall's 
correlation and partial-correlation coefficients (if available), along with significance levels for these 
three samples. We present the statistical results of Sample A for comparing the correlation results 
of C IV and Fe Ka. Because the combined sample contains censored data for the a ox values, we 
must use survival analysis to calculate the correlation coefficients. However, algorithms are not 
available for calculating all the correlation and partial-correlation coefficients for censored data. 
For example, the empty entries in Table[5]are due to the unavailability of corresponding algorithms. 
In order to compare the changes of correlation strength when the third parameter is controlled, we 
construct a reduced sample with the censored data removed, considering that this only excludes 
a small fraction (~ 6%) of the entire data set and should not affect the statistical properties of 
the sample. We can see that the C IV BEff is significantly weakened when a ox is held fixed; the 
correlation coefficient drops from —0.580 to —0.224 (Spearman). On the other hand, the correlation 
coefficient between EW(C iv) and a ox also drops significantly when Zj,(2500 A) is held fixed, from 
0.615 to 0.332. This implies that both a ox and Z i ,(2500 A) are driving the change of EW(C iv). 

The Fe Ka BEff plot (not show n) of our sample (Sample A) is very similar to the correlation 



shown of Fig. 4 in lJiang et al.l (|2006l ) except that our sample size is smaller. Fig. [5]shows EW(Fe Ka) 
plotted against a ox . The Spearman test gives a much weaker correlation coefficient (—0.230 with 
Po = 0.111) than that for C iv (—0.304 with P < 0.001). In addition, simple \ 2 fitting produces a 
slope of 0.046 ± 0.154, consistent with zero. The consistency between the correlation analysis and 
regression result provides strong evidence that EW(Fe Ka) is not correlated with a ox . 



4.3. Effects of AGN Variability on BEff Relation Scatter 

Because of the ubiquity of AGN variability, combined with the different times of the optical 
and X-ray observations, our values of a ox do not reflect the spectral hardness at a specific time 



} Po is the confidence level of the null hypothesis. Therefore, the smaller Po is the more likely the correlation exists. 
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but are randomly distributed around their mean values. The deviation of a ox from its mean value 
would be ~ 0.083, assuming that the variation amplitudes are 30% for ^(2500 A) and 40% for 
l v (2 keV) (e.g., Strateva et al. 2005, GBS08). 

To check how much of the scatter of our correlations could be attributed to variability, we 
performed two simple tests on our combined sample. We follow the method used in § 3.1 of GBS08 
and introduce AlogEW, which is the difference between the observed EW(C iv) and the EW 
calculated from linear regression (Eq. ED, i.e., AlogEW = log EW - log EW(Z„ (2500 A)). We 
define \i and a as the mean and disp ersion of the distributio n of A log EW. To calculate \x and a, 
we maximize the likelihood function (jMaccacaro et al.lll988l ): 



L 



n 



i 



27T (af + a") 



: exp 



(AlogEW,; 



(6) 



in which the subscript i represents each object and <7, is the uncertainty of AlogEW associated 
with each AlogEW,. The maximization of L requires fi = 0.01 and a = 0.23. 

Next, we estimate the potential scatter of AlogEW due to variability. We need to consider 
two terms: logEW and log EW(/„(2500 A)). To first approximation, EW ~ /iine//cont in which 
/line is the emission-line flux and the / CO nt is the continuum flu x. The emis s ion-lin e variability 
of six luminous quasars at z = 2.2-3.2 was recently reported by iKaspi et al.l (|2007l ). The mean 
fractional variation K^J^I is ~ 0.096 by averaging the fractional variation of C IV A1549 of all 
six quasars. The C IV emission-line variability of a number of Se yfert galaxies has been studied, 

NGC 7469 



i ncluding Fairall 9 ( Rodritmez-Pascual et all Il997fl. N GC 5548 (jClavel 



et al 



1991 



(jWanders et al.lll997f ). NGC 3783 (jEeichert et al.lll994l ). and 3C 390.3 (IQ'Brien et al.lll998h . By 
averaging the fractional variations of Seyferts and quasars above, we obtain an average emission- 
line variation (-F var ) = 0.130. The typical variation of /^(2500 A) is ~ 30% (e.g., Strateva et al. 
2005). Therefore, the scattering of logEW contributed from variability is estimated (assuming all 
independent variables are Gaussian distributed) as 



a (A logEW) 




0.144. 



In the calculation above, a = 0.198 (Eq. ([7])). This exercise indicates that at least 
scatter around the BEff in our sample can be attributed to AGN variability. 



of the 



We performed a similar test for the a ox -EW(C iv) relation. Because the set of a ox contains 
censored data, we can only use the reduced data set (258 objects). The maximization of L (Eq. [6]) 
yields fi = 0.01 and a = 0.22. The potential dispersion of this relation assuming all scatter comes 



is defined as the RMS of the intrinsic variability relative to the mean flux IjRodriguez-Pascual et al.lll997f ) 
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from variability is estimated (assuming Gaussian distributions) as 
a(AlogEW) 



1 



In 10 



line 



/line 



+ 



fa 



+ (0.3838a) 2 



<H„(2500 A) 



^(2500 A) 



+ 



5l v {2 keV) \ 2 ' 
l v {2 keV)J 



In the calculation above, a = 1.035 (Eq. [5]). This indicates that variability produces at least 75% 
of the scatter around the a ox -EW(C iv) relationship. 

In summary, the above two tests demonstrate that a substantial fraction, if not the majority, 
of the scatter in the correlations above can be attributed to X-ray and UV/optical variability. It 
should be possible to make the correlations tighter if the UV/optical and X-ray data are observed 
simultaneously. 



0.163. 



4.4. Regressions of EW and Luminosity 

The BEff provides a potential avenue to infer the luminosity of a quasar from emission-line 
observations. Type la supernovae (SNe) are treated as classical standard candles (e.g., Phillips 
1993; Burrows 2000), but only a few are observed beyond z ~ 1.5. If quasars, which are much 
easier to detect than SNe and can be observed to much a higher redshift, can be used as standard 
candles, they would be an important tool for cosmological studies. Soon after the discovery of the 
BEff, many investigations considered the possibility of treating emission-line EW as a luminosity 
indicato r (e.g., Baldwin 1977b: Wampler 198 0). Unfortunately, the C iv BEff usually has a large 



scatter (jOsmer fc Shields! 11999k IShieldsl 120071 ) ; given the small slope in the logEW-logZ„ plot (on 
the order of —0.2), the predicted luminosities are very inaccurate. It has also been shown that the 
C iv BEff is redshift dependent (Francis & Koratkar 1995; Vanden Berk et al. 2009, in preparation), 
making it a less valuable probe of cosmology. 

Because we are focusing on the influence of a ox at the moment, and will put the redshift factor 
aside, we will concentrate on the issue of whether the scatter can be reduced if we regress EW(C iv) 
with ^(2500 A) and/or a ox , and if the prediction of luminosity can be made more accurate with 
this approach. The linear-regression results of EW(C iv) with a ox are already shown in Eq.([5]). 
Similar regressions from ^(2500 A) and both of Z„(2500 A) and a ox , using the fully parametric EM 
algorithm, are 

log EW(C iv) = (-0.198 ± 0.015) log /^(2500 A) + (7.764 ± 0.461) (7) 
logEW(C iv) = (-0.107 ±0.021) log ^(2500 A) + (0.615 ± 0.106)a ox + (5.944 ± 0.536). (8) 

in which EW(C iv) is in A and /^(2500 A) is m erg s _1 Hz" 1 . The last regression was performed 
on the combined sample without the censored data (258 objects) because the EM algorithm in 
ASURV does not allow both independent variables to contain censored data points. 



To evaluate the scatter, we subtract the predicted EW values calculated using the above 
equations from the observed values and compute the RMS values of the residuals. We see a slight 
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improvement of the RMS values, using a ox and both / l/ (2500 A)+q ox (Table [6|). The last regression 
result (Eq.®) is consistent with Eq.(ll) in GBS08, indicating that at least part of the scatter of 
the BEff is due to a ox . 

Next, we regress luminosity against EW(C iv) and/or a ox , using the combined sample. The 
results are 

log ^(2500 A) = (-1.980 ± 0.154) log EW(C iv) + (34.076 ± 0.256) (9) 
log Z„(2500 A) = (-0.852 ± 0.168) log EW(C iv) - (2.864 ± 0.262)a ox + (27.668 ± 0.635). (10) 

The last regression was performed on the censored-data-excluded sample. We then calculate the 
RMS values of the residuals after subtracting the predictions from the equations above (Table [6]). 
The RMS value shrinks by 18% using EW(C iv)+a ox , compared to using EW(C iv) alone. We 
use the standard -F-test to check if the two sets of residuals have consistent variance. The testing 
gives an F-statistic of 1.48 with a significance 0.002, indicating that these two sets of residuals have 
significantly different variances and 18% is a statistically significant improvement. To use quasars 
as standard candles via the BEff, we should at least confine the luminosity within an uncertainty of 
30%, or, equivalently, rms < 0.1. This cannot be achieved using our current dataset and controlled 
parameters. 



5. Relation Between Fe Ka and C iv 

Fe Ka is important in AGN studies because it is the strongest emission line appearing in 
the X-ray band. However, the strength of this emission line varies significantly from object to 
object, and the line is not detected in most X-ray observations of quasars. Given this lack of direct 
observational measurement, it would be useful to develop a way to predict the expected strength 
of the line empirically. 

The EW(C iv) and EW(Fe Ka) do not exhibit a significant correlation (Fig. [6]) in Sample A; 
the correlation has a low Spearman rank-correlation coefficient (0.319 with Pq = 0.027). Although 
the EWs of the two lines are not well correlated, their luminosities and fluxes are strongly corre- 
lated (Fig. [7J), with Spearman correlation coefficients 0.529 (P < 0.001) and 0.551 (P < 0.001), 
respectively (Table [7]). The linear correlations regressed in Fig. [7] are 

log L(Fe Ka) = (0.588 ± 0.079) log L(C iv) + (16.164 ± 3.416) (11) 
log /(Fe Ka) = (0.978 ± 0.188) log /(C iv) - (2.082 ± 2.228). (12) 

One must always question the significance of relations such as (jllh because even if there is no 
correlation between the observed fluxes of the lines, the fact that the line luminosities for a given 
object contain the same distance factor will introduce an apparent correlation in the luminosities. 
To investigate whether this effect is important for our study, we perform a test in which we conduct 
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correlation and regression analysis for a sub-sample of Sample A. Objects in this sub-sample have 
similar redshifts, and thus they have approximately the same distances. First, we use objects with 
0.06 < z < 0.09 because this redshift bin contains a large number of objects. This sub-sample 
contains 10 objects. The correlation coefficient of /(Fe Ka)-/(C iv) is 0.624 (Pq = 0.061) and of 
L(Fe Ka)-L(C iv) is 0.709 (P = 0.033). The regression results are 

log L(Fe Ka) = (0.748 ± 0.200) log L(C iv) + (9.330 ± 8.702) (13) 
log /(Fe Ka) = (0.771 ± 0.241) log /(C iv) - (4.280 ± 2.819). (14) 

Both the correlation and regression results of this sub-sample are consistent with the results of 
the entire sample within uncertainties, suggesting that the luminosity correlation of the C iv and 
Fe Ka lines may be real and is not a consequence of multiplying the two fluxes of a given object 
by the same large distance factor. The luminosity of C iv emission line increases faster than the 
luminosity of Fe Ka. 

One might be concerned that the correlation between L(Fe Ka) and L(C iv) is artificial 
because the calculation of the Fe Ka measurements involves L(2-10 keV) (Eq. (pi])), which is 
proportional to l u {2 keV) (Eq. (jTJ) ) , and l v {2 keV) is correlated with ^(2500 A), which is correlated 
with EW(C iv), i.e., the C iv BEff. EW(C iv) is calculated from continuum and emission line 
luminosity, so apparently, the L(Fe Ka) and L(C iv) are not independent before we perform the 
co rrelation. However , our calculation of Fe Ka is simply reversing the .F(2-10 keV) /EW calculation 



of ljiang et al.1 (j2006T ). so F(Fe Ka) is not actually dependent of F(2-10 keV) and hence Z„(2500 A). 
In essence, we have the values of F(Fe Ka) independent of jF(C iv). Therefore, our L(Fe Ka) and 
L(C iv) correlation, which is expected from existing relations, is not an artifact correlation. 



That the EWs of C IV and Fe Ka are uncorrelated is consistent with the result by lPage et al 



((2004J) and further demonstrates that the C iv and Fe Ka emission lines are unlikely to have 
the same origin. This result is not surprising because Fe Ka and C IV are produced in different 
processes. The correlation between their luminosities is probably a combination of effects between 
their EWs (uncorrelated) and continuum (strongly correlated). The flux correlation, although 
empirical and not very tight, is a useful first order estimation of the Fe Ka line flux given the UV 
spectra in the rest-frame of an AGN. 



6. Discussions and Conclusions 

We have compiled a sample of 272 Type 1 AGNs and quasars that have UV and X-ray mea- 
surements, among which Fe Ka emission lines are detected in 50 objects. The sample covers a wide 
range of redshift (0.009 < z < 4.720), and a wide range of luminosity from Seyfert galaxies to the 
most-luminous quasars (27.81 < log Z„ (2500 A) < 33.04). These properties allow us to study the 
overall properties of AGNs rather than focusing on a particular redshift or luminosity. It also has 
a high X-ray detection rate (~ 96%), which lets us obtain robust statistics. We have performed 
correlation and regression analyses using this sample and draw the following conclusions: 
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1. The C iv BEff is driven by both a ox and ^(2500 A), or equivalently, by l u (2 keV) and 
Zj,(2500 A). This implies that changes in the ionizing flux induce changes in the ionization 
state of the BELR, producing more C iv ions when the SED becomes harder and vice versa. 
This is supported both by correlation and regression anlayses: 

• The partial correlation between EW(C iv) and ^(2500 A) when a ox is controlled is 
weaker than the regular correlation between EW(C iv) and /^(2500 A). 

• The scatter in the linear regression decreases when we regress EW with a ox +l u {2500 A) 
compared with ^(2500 A) alone. 

Although the reduction of the scatter due to adding another regression parameter is not 
sufficiently large to treat quasars as standard candles, it demonstrates that a significant 
fraction of the scatter attributes to a ox , and can be reduced by including it in regression 
analysis. 

2. EW(Fe Ka) exhibits no strong correlation with either a ox or EW(C iv). This implies that 
Fe Ka is not likely to have the same origin as C iv. 

3. There may be a correlation between the luminosities of Fe Ka and C iv with a logarithmic 
slope of 0.588 ±0.079. This correlation is possibly because both of these two quantities involve 
a factor related to the scale of the line emitting regions and the slope indicates that the C iv 
emission line luminosity increases faster than the Fe Ka. 



Although a ox is a fundamental influence on EW(C iv), there is still a significant scatter in the 
EW(C IV) -a ox diagram. As we have demonstrated, most of the scatter is contributed by variability, 
but another likely contribution source is the nature of a ox which only connects the flux points at 
2500 A and 2 keV but misses the Big Blue Bump, which is expected to play an important role 
in the photoionization process. The shape of an AGN SED ca n be very different depend ing on 
Eddington ratio (-Z/boi/^Edd) but still have a fairly constant a ox (jVasudevan fc Fabianl 120071 ). It is 
perhaps more appropr iate to use a point near ~ 250 A instead of 2500 A to calculate a revised a ox 
(jShemmer et al.ll2008l ). This new a ox might be more strongly correlated with EW(C iv). However, 
this requires challenging observations that cannot be achieved at present for most AGNs. 
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Fig. 1. — The luminosity-redshift diagram of the combined sample containing 272 objects distin- 
guished by observational facilities. The gap for 0.5 < z < 1.5 arises because of the detection limit 
of HST/IUE and wavelength coverage of the SDSS. The dashed line marks the position where 
Z„(2500 A)= 30.5. 
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Fig. 2. — Continuum and C iv emission-line fit examples. In each panel, the upper spectrum 
is the original and the lower spectrum is continuum subtracted; blue solid curves are the fits 
to the spectra. The cyan dotted curves are emission-line components. For PG 0947+396 and 
3C 120, these components include the powerdaw continua and two Gaussian profiles for C iv. 
For SDSS J100129. 64+545438.0, we also plot the iron emission forest and small Balmer bump 
components, although they are so weak that they are almost invisible. The rest-frame spectral 
resolutions of these spectra are, from top to bottom, ~ 1.6 A, ~ 5 A, ~ 0.8 A. 
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Fig. 3. — The C iv BEff diagram for the combined sample. The EW(C iv) is given in A, and 
/^(2500 A) is in units of erg s 1 Hz 1 . Points are distinguished by luminosity; the short solid line 
is the best linear fit to low-luminosity points and the long one is the best linear fit to the entire 
sample. 
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Fig. 4. — The correlation between EW(C iv) and a ox . Upper limits on a ox are marked with arrows. 
The solid line is the best linear fit using the EM algorithm. 
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Fig. 5. — The correlation between EW(Fe Ka) and a ox for the core sample. The solid line is a 
linear fit using the ASURV package for the censored data. The typical error bar of the data is 
displayed at the top-right corner of the plot. 
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Fig. 6.— Plot of log[EW(Fe Ka)/E (Fe Ka)] vs. log [EW(C iv)/A (C iv)] for the core sample. 
Because the units of EW(Fe Ka) and EW(C iv) are different, we divide them by central energy 
Eo(Fe Ka) = 6.4 keV and central wavelength Aq(C iv) to make them dimensionless. 
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Fig. 7. — Plot of logL(Fe Kq) vs. logL(C iv) (upper panel) and log/(Fe Ka) vs. log/(C iv) 
(lower panel) for the core sample. Upper limits are denoted as downward arrows. The solid lines 
are the best linear fits to the data using ASURV. For comparison, we show the best linear fits with 
a unity slope in dotted lines. The data symbols and typical error bars are labelled at the bottom 
right corners. 
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Table 1. Summary of samples. 





Sample 


SDSS 


HST 


IUE Total 


Rcdshift range 


log i„(2500 A) range 






A 







34 


16 50 


0.009-1.735 




27.81-31.69 






B 




98 





98 


1.7-2.7 




30.53-31.67 






C 




91 


13 


20 124 


0.015-4.720 




28.12-33.04 






Combined 


189 


47 


36 272 


0.009-4.720 




27.81-33.04 








Table 2. 


UV properties of 


the combined sample. 




Object 




2 


logi„ 


(2500 A) 


EW(C IV) 




Flag 1 


UV/Optical 


Sample 








(erg s 


- 1 Hz" 1 ) 


(A) 






Instrument 




NGC4593 




0.009 


27.81 ± 0.15 


141.0 ± 23.2 


-1.008 ± 0.066 





HST 


A 


NGC3783 




0.010 


28.40 ± 0.08 


144.0 ± 9.4 


-1.251 ± 0.045 





HST 


A 


Mkn352 




0.015 


28.12 ± 0.15 


200.0 ± 12.9 


-1.326 ± 0.067 





IUE 


C 


Mrkl044 




0.016 


28.63 ± 0.04 


60.1 ± 1.3 


-1.523 ± 0.037 





HST 


A 


NGC7469 




0.016 


28.72 ± 0.28 


129.0 ± 25.1 


-1.319 ± 0.111 





HST 


A 


MCG8-11- 


-11 


0.020 


28.34 ± 0.19 


233.0 ± 39.1 


-1.182 ± 0.082 





IUE 


C 


Mkn79 




0.022 


28.75 ± 0.09 


135.0 ± 11.3 


-1.239 ± 0.049 





IUE 


C 


Mrk335 




0.026 


29.29 ± 0.18 


75.5 ± 7.7 


-1.503 ± 0.077 





HST 


A 



Note. — Tablc[2]is published in its entirety in the electronic version of the Astrophysical Journal. The portion is 
shown here for guidance regarding its form and content. 

"""A value of "1" indicates that a ox for this object is an upper limit. 



Table 3. UV and X-ray properties of Sample A objects. 



Object 


z 


EW(Fe Ka) 
(eV) 


logi„(2 kcV) 
(erg s _1 Hz -1 ) 


log L(Fe Ko 
(erg s _1 ) 


) log/(FcK Q ) 
(erg s — 1 cm - 2 ) 


log L (Civ) 
(erg s _1 ) 


log /(Civ) 
(erg s — 1 cm 2 ) 


NGC4593 

NGC3783 

NGC7469 

Mrkl044 

Mrk335 


0.009 
0.010 
0.016 
0.016 
0.026 


qs n+ 210 
120 n+ 14 

105.0+.=!;° 

< 54.0 


25.18 ± 0.09 
25.14 ± 0.09 
25.28 ± 0.09 
24.66 ± 0.09 
25.38 ± 0.09 


41 05 + °" 09 

4i.i8±° ;° l2 

40.81+g-" 
< 40.99 


, O1 +0.08 

— ■"•■"■-n.ii 
-12.981°;- 

< -13.19 

,-, sd +0.22 
-11.84_ 4g 

< -14.42 


41.70 ± 0.03 
42.50 ± 0.01 
42.35 ± 0.03 
42.02 ± 0.01 
42.75 ± 0.02 


-11.55 ± 0.03 
-10.82 ± 0.01 
-11.43 ± 0.03 
-11.76 ± 0.01 
-11.43 ± 0.02 


Mrk590 
Mrk290 
Mrk493 


0.026 
0.030 
0.031 


121 n+ 650 

eo c + 56.5 
t>s - t >_43.4 
< 101.0 


25.25 ± 0.09 
25.62 ± 0.09 
25.31 ± 0.09 


41 21+ 019 

41 27+ 029 
* x "" -0.59 

< 41.19 


-~3 41+°- 17 
"•"-0.28 

-12 62+ 08 
< -13.55 


42.56 ± 0.01 
42.79 ± 0.01 
42.21 ± 0.03 


-11.64 ± 0.01 
-11.52 ± 0.01 
-12.15 ± 0.03 


Note. — 


Table [3] is 


published in its entirety in the electronic version 


of the Astrophysical Jot 


irnal. The port: 


ion is shown here for 



guidance regarding its form and content. 



Table 4. Hypothesis and linear fitting results. 



Fig 


Sample 


X 


y 


p 1 (^o 2 ) 


k 3 


m 


Combined 


log ! A (2500 A) 


logEW A (C IV) 


-0.559(< 0.001) 


-0.198 ± 0.015 


m 


Combined with log i„ (2500 A) < 30.5 


log l x (2500 A) 


logEW A (C IV) 


-0.465(< 0.001) 


-0.218 ± 0.048 


m 


Combined 


a ox 


logEW A (C IV) 


0.755(< 0.001) 


1.035 ± 0.075 


m 


Sample A 


C* x 


logEW eV (Fe Ka) 


0.104(0.470) 


0.247 ± 0.230 4 



Spearman rank correlation coefficient. 

2 Significance levels of Spearman's rank correlation. 

3 Slope from EM algorithm. 

4 C alculated for Sample A using the Buckley- James method in the ASURV software package ( Laval ley, Isobc &i Fcigclson 
119921 1. 
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Table 5. Correlation and partial-correlation analysis results. 



X 


y 


2 a Spearman 


Pearson 


Kendall 


Combined sample (279) 


logEW(C IV) 
logEW(C IV) 


log It, (2500 A) 
log !„ (2500 A) 


a ox ■ ' ■ 
-0.599(<0.001) 






-0.237 b 
-0.400(<0.001) 


a os 


logEW(C iv) 
logEW(C IV) 


log l„ (2500 A) 

0.607(<0.001) 






0.260 






Comb 


incd Sample Without Censored Data 


(258) 






logEW(C IV) 
logEW(C IV) 


log i„ (2500 A) 
log l u (2500 A) 


a ox -0.224 
-0.580(< 0.001) 




-0.311 
-0.605 


-0.220 
-0.417(< 0.001) 


a ox 


log EW(C IV) 
log EW(C IV) 


log l v (2500 A) 0.332 
0.615(< 0.001) 




0.258 
0.587 


0.284 
0.450 




Sample A (49) 


logEW(C IV) 
logEW(C IV) 


log i„ (2500 A) 
log l v (2500 A) 


a ox 0.120( 0.291) 
-0.304(0.001) 


0.074( 
-0.304( 


0.214) 
0.001) 


0.036 
-0.180( 0.001) 


logEW(Fe Ka) 
logEW(Fo Ko) 


log(„(2 kcV) 
log l„{2 kcV) 


a ox .-. c 
-0.230( 0.111) 






< 0.001 


a ox 

Otox 


logEW(Fe Ka) 
log EW(Fe Ka) 


log ;„(2 kcV) 

0.104( 0.470) 






-0.130 



Note. — The number in the parentheses after a correlation coefficient is its significance level. Because we are 
considering the null hypothesis, a small number indicates a possibility of a strong correlation. 

a z is the controlled parameter. If a 2 entry is not empty, we calculate the partial correlation of x and y while 
controlling for 2; otherwise, we only calculate the correlation of x and y. 

^For the Kendall's partial r correlation coefficients, generally, the sampling distribution is unknown; therefore, 
the probability values are not available (Kendall 1938, 1970). 

c,c ^Because EW(Fe Kck) is a set of censored data, we have to use a generalized correlation statistics that can 
deal with censored data to derive the coefficients. For the non-partial correlations, we used the ASURV software 
package (Lavallcy, Isobc Sz Fcigclson 1992) which only provides Spearman's p and Kendall's r correlation correlation 
coefficients. For the PCA, we only use the Kendall's partial r correlation statistics jAkritas &c Siebeit|fl996| ) . 



Table 6. The RMS values of residuals after regression from different variables. 



Independent variables 


Dependent variable 


RMS values 


Z„(2500 A) 


EW(C iv) 


0.231 


CTox 


EW(C iv) 


0.228 


«„(2500 A)+a ox 


EW(C iv) 


0.217 


EW(C iv) 


Z„(2500 A) 


0.747 


CTox 


Z„(2500 A) 


0.645 


EW(C iv)+a ox 


Z„(2500 A) 


0.615 



Note. — To consistently compare the RMS values, we com- 
pute them using the combined sample without the censored 
data. 
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Table 7. Correlation and regression analysis for EW, emission line luminosity, and flux data 

between Fe Ka and C iv. 



Relations 


P(-P) 1 


k 2 


EW 


0.319( 0.027) 


0.291±0.131 


L(line) 


0.529(<0.001) 


0.588±0.079 


/(line) 


0.551(<0.001) 


0.978±0.188 



1 Correlations are tested using Spearman's 
p and the significance level (P) is evaluated 
against the null hypothesis. 

2 We use the Buckley-James method to do 
linear regression and k is the slope. The com- 
putation is done using the ASURV softw are 
package jLavallev, Isobe fe Feiee lson 1992) 



